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■ Abstract 

We present the mechanism of the dynamical supersymmetry breaking at the metastable vacuum 
recently uncovered in the J\f — 1 U(N) supersymmetric gauge theory that contains adjoint 
superfields and that is specified by Kahler and non-canonical gauge kinetic functions and a 
superpotential whose tree vacua preserve M = 1 supersymmetry. The overall U{1) serves as 
the hidden sector and no messenger superfield is required. The dynamical supersymmetry 
breaking is triggered by the non-vanishing D term coupled to the observable sector, and is 
realized by the self-consistent Hartree-Fock approximation of the NJL type while it eventually 
brings us the non-vanishing F term as well. It is shown that theoretical analysis is resolved 
as a variational problem of the effective potential for three kinds of background fields, namely, 
the complex scalar, and the two order parameters D and F of supersymmetry, the last one 
being treated perturbatively. We determine the stationary point and numerically check the 
consistency of such treatment as well as the local stability of the scalar potential. The coupling 
to the M = 1 supergravity is given and the gravitino mass formula is derived. 
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1 Introduction 



Spontaneous breaking of rigid supersymmetry occurs much less frequent compared with that of 
internal symmetry in quantum field theory and has attracted much interest [HE] of theorists for 
over the three decades. Mass hierarchy in elementary particle physics indicates that it is most 
desirable to break M = 1 supersymmetry dynamically. In fact, under the non-renormalization 
theorem [3], no holomorphic operator is generated in perturbation theory and instanton gener- 
ated nonperturbative superpotentials have been the major source of dynamical supersymmetry 
breaking (DSB). 

In this paper, we focus our attention on general rigid M = 1 theory in four spacetime 
dimensions consisting of vector superfields and chiral superfields in the adjoint representation 
which permits a non-canonical gauge kinetic function r a b (that may follow from the second 
derivative of the prepotential) and hence the D term-gaugino-matter fermion (or D term-Dirac 
gauginos) nonrenormalizable coupling. It has recently been shown in refs. [HE] that, in this 
general situation, supersymmetry is dynamically broken in the metastable vacuum. The mech- 
anism that triggers the DSB is the condensate of the Dirac bilinear above, forcing one of the 
order parameters D of supersymmetry to be non-vanishing. This is very much reminiscent of 
the Nambu-Jona Lasinio (NJL) theory [6j[7] of broken chiral symmetry and hence the BCS 
superconductivity [H[9], being formulated in terms of the effective action of the auxiliary field 
whose stationary value is the order parameter. The method of approximation employed is the 
self-consistent Hartree-Fock approximation where the tree and the one-loop contributions are 
regarded as comparable. Once this mechanism operates, non-vanishing F term is shown to be 
induced and contributes, for instance, to the mass of the fermions. The mechanism requires 
massive adjoint scalars, in particular, the scalar gluons and, together with the feature that the 
D term triggers the breaking, is quite distinct from the previous proposals [TOHM] of DSB both 
from theoretical and experimental perspectives. The overall U(l) where the non-vanishing D 
and the Nambu-Goldstone fermion (NGF) reside serves as the hidden sector and no messenger 
field is necessary [3] as non- vanishing third prepotential derivatives connect the U(l) sector 
with the observable SU(N) sector [T5l[T6] . 

While our treatment of the theory bears much resemblance with that of the NJL theory, 
there is one important complexity which has no counterpart in the NJL and which we did 
not emphasize in ref. pQ. In NJL, aside from the pseudoscalar auxiliary massless singlet field 
commonly denoted by ir, there is only one singlet auxiliary scalar field denoted by a in the 
effective action, which is the order parameter of chiral symmetry. (See appendix A.) In other 
words, the stationary condition of energy with respect to the scalar is at the same time the 
stationarity with respect to the order parameter (the gap equation). This is not the case here. 
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After treating the U(N) singlet real auxiliary field F as perturbation, we have two kinds of 
background fields in the effective potential: these are the singlet complex scalar if and the 
singlet auxiliary field D. The stationarity of energy with respect to the scalar and that with 
respect to the order parameter are one and the other and both must be imposed simultaneously. 
In this paper, we will mainly deal with such multi-variable variational problem in depth and 
present the solution which is the local minimum of the scalar potential. We will also include a 
few other materials which have phenomenological implications. We work in the unbroken phase 
of U(N) and invoke U(N) invariance of the expectation values to suppress indices often. 

In the next section, we start out from exhibiting the component action from that of the 
superspace, state the set of assumptions we have made in [HE] and in this paper and give the 
Noether current associated with M = 1 rigid supersymmetry. We review the original reasoning 
that has led us to the D-term triggered dynamical supersymmetry breaking. We set up the 
background field formalism to be used in the subsequent sections, separating the three kinds 
of background from the fluctuations. The action can be coupled to M = 1 supergravity and 
we derive the gravitino mass formula via the super-Higgs mechanism associated with the non- 
vanishing D-term. The action contains a sequence of special cases in which the gauge coupling 
function and the superpotential are related in a specific form, including the one where the rigid 
M = 2 supersymmetry is partially broken to M = 1 at the tree level [T5|[T7]. In section three, we 
elaborate upon our treatment of the effective potential with the three kinds of background fields 
as well as the point of the Hartree-Fock approximation in refs. (HE]. Section four is the main 
thrust of this paper. We present our variational analyses of the effective potential in full detail. 
Treating one of the order parameters F as an induced perturbation, we demonstrate that the 
stationary values (£)*, tp*, tp*) are determined by the intersection of the two real curves, namely, 
the simultaneous solution to the gap equation and the equation of <p stationarity (the energy 
condition). Numerical analysis is provided that demonstrates the existence of such solution 
as well as the self-consistency of our analysis. The second variation of the scalar potential is 
computed and the local stability of the vacuum is shown from the numerical data. We finish 
our paper with summary and brief comments on the issue of regularization and subtraction 
schemes. 

In two of the appendix on rudimentary materials to be referred to in the text, we take a brief 
look at the NJL effective action and recall a formula of the second variation of a multivariable 
function. Phenomenological applications of our finding and the estimate of the longevity of our 
metastable vacuum have been given in [HE], which we do not repeat in this paper. 
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2 The action, assumptions and some properties 



The action we work with in this paper is the general M = 1 supersymmetric action consisting 
of chiral superfield $ a in the adjoint representation and the vector superfield V a with three 
input functions, the Kahler potential K(Q a , $ a ) with its gauging, the gauge kinetic superfield 
T"afe($ a ) that follow from the second derivatives of a generic holomorphic function J-"($ a ), and 
the superpotential W(<& a ). 

C = j d 4 9K($ a ,$ a ) + (gauging) + J d 2 9lm^r ab (^ a )W aa W b a + (J d 2 6W($ a ) + c.c.^(2.1) 

The gauge group is taken to be U (N) and, for simplicity, we assume that the theory is in the 
unbroken phase of the entire gauge group, which can be accomplished by tuning the superpo- 
tential. We also assume that third derivatives of J-"($ a ) at the scalar vev's are non-vanishing. 

2.1 action and component expansion 

The component Lagrangian of eq. (12.11) reads 

Kahler £gauge ^sup; (2-2) 

where 

^Kahler = g ab V^ a V»f - '-g^W^ + % -g ah V'^ a o^ + g ab F a F b 

~g ab - c F a ^ c - X -g hc ^F c r^ + -^g ab (\ c rK b + WK a ) + \D a D a , (2.3) 

£gauge = -^ AV V 6 - l^abD^Jf ~ ^) ab F^F^ - 1 (^) ab e^ F^ 

\T ahc ro v o^\ a - F abc \ a a^ c )F b 



+ l -(^) ab D a D b + ^-(T abc ^\ a + T ab ^ c \ a )D b + l -F abc F c \ a \ b - l -F abc F c \ a \ b 
- l -J" abcd ^ d \ a \ b + l -T ab cdi> c i> d \ a \ b , (2.4) 

o o 

£su P = F a d a W - l -d a d b Wr^ h + c.c, (2.5) 



where 

Va = -\(F b f b a & + F b ft<P C ) (2.6) 

and f b c is the structure constant of SU(N). Note that an equation of motion for F a is F a = 
—g ab d b W + fermions. We also assume (F a ) trcc = —(g ab d b W) tICC = at the tree level. At the 
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lowest order in perturbation theory, there is no source which gives vev to the auxiliary field D°: 
(-D°)tree = 0. The U(N) gaugino is massless at the tree level while the fermionic partner of the 
scalar gluon receives the tree level mass m a = m = (g 00 dod W) t ree- 

2.2 assumptions 

While we have already stated, it is useful to recapitulate here a set of assumptions made in order 
to address better the question of dynamical supersymmetry breaking within our framework. 

1) a general M = 1 supersymmetric action of chiral superfield $ a in the adjoint representa- 
tion and the vector superfield V a with the three input functions, namely, the Kahler potential 
K ($ a , <5 a ) with its gauging, the gauge kinetic superfield r ab ($ a ) that follow from the second 
derivatives of a generic holomorphic function J-"($ a ), and the superpotential W(Q a ). 

2) third derivatives of J r ($ a ) at the scalar vev's are non-vanishing. 

3) the superpotential at tree level preserves M = 1 supersymmetry. 

4) the gauge group is U(N) and the vacuum is taken to be in the unbroken phase of U(N). 
It is in principle straightforward to extend this to the (partially) broken cases where U(N) is 
broken into the product groups. The variational analyses we carry out in section four, however, 
become more complex and we will not address this in this paper, See the comment at eq. A3. 22[) 



2.3 supercurrent 

We give here an off-shell form of the Af = 1 supercurrent. 

mSM" = V2g ab r )l <7 v a^ a Vj b + V2ig abm a^ a F b 

-iT abVl a u \ a F^ b + ^F ab e^o 5 r ]1 a v \ a FP Sh - % -F ab r ]1 a»\ a D b 

+^ (T abc r^^X b - T ab ;\ c o»o^ h ) T)\a v \ a . (2.7) 
Equations of motion for auxiliary fields are 

D a = ~g ab 2 b - ^=g ab (J W ^A C + T bcd ^p d X c ) , 

F a = -g^dtW - % -g ab {T bcd ^ d ~ J r b cA c \ d ) ■ (2.8) 

Once the U (N) invariant components of the auxiliary fields, D° and F° receive non- vanishing 
vev's together with U (N) invariant scalar vev's, the second and the fifth terms of the RHS of 
eq. (12. 7p at these vev's clearly develop a one-body fermionic operator non-vanishing at zero 
momentum [T8H20] : this particular combination of tp° and A creates the one-particle state 
which is identified with the Nambu-Goldstone fermion [T9l. 
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2.4 original reasoning of DDSB 

In ref. jl], it was shown that the vacuum state of the theory, albeit being metastable, develops 
a non-vanishing vev of an auxiliary field D° in the Hartree-Fock approximation. The theory, 
therefore, realizes the D-term dynamical supersymmetry breaking. The relatively simple es- 
timate has shown that the vacuum can be made long lived. Let us recall a few more key 
aspects. 

The part of the lagrangian which produces the fermion mass matrix of size 2iV is 

-W)( ° -^ Db )( X °) + (c,.). (2.9) 

2 '\-^t^d* d a d,w 

It was observed that the auxiliary D a field, which is an order parameter of M = 1 super- 
symmetry, couples to the fermionic (but not bosonic) bilinears through the third prepotential 
derivatives: the non-vanishing vev of D° immediately gives a Dirac mass to the fermions. Eq. 
(EHD implies 

(D°) = -7^(9°° (-WA C + Jo ca #A c )>, (2.10) 

telling us that the condensation of the Dirac bilinear is responsible for (D°) ^ 0. 
We diagonalize the holomorphic part of the mass matrix: 

{-^(FoaaD ) (8 a d a W) )' ^ ' } 

Note that the non-vanishing third prepotential derivatives are J-oaa where a refers to the gen- 
erators of the unbroken gauge group. By an orthogonal transformation, we obtain the two 
eigenvalues of eq. (12. lip for each generator, which are mixed Majorana-Dirac type : 



Introducing 



we obtain 



^^whW'^i' (2 - 12) 



|A£i| 2 = |<a„W||A<fJ| 2 . (2.14) 



It was also shown in ref. [1] that the non-vanishing F° term is induced by the consistency 
of our procedure of computation. (See also [2T1I27]). This is because the stationary value of the 
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scalar fields gets shifted upon the variation (the vacuum condition). The final mass formula for 
the SU (N) fermions is to be read off from 

= -\(g a,a)(F°)rr + \{ ^aa) (F°> A a A" - \ (d a d a W) ^ + ^ ( ^Oaa) P ^ (D°) 

1 ^ f\ a (r)\ 
= - 5 E *(x) at Ma, a y a (x), * a (x) = [ y 1 . (2.15) 

We will write down the explicit form in the next subsection. See eqs. (I2.17p . (I2.18p . ( 12 . 1 9[) and 
(12.201) . A main remaining point is how to establish the procedure in which the stationary values 
of the scalar fields, D° and F° perturbatively induced are determined, which we will resolve in 
this paper. 



2.5 quadratic part of the quantum action 

In this subsection, we write down parts of the action with the background fields for the com- 
putation of the one-loop determinant in the next section. The linear terms that arise upon 
separation into quantum fields and background fields are dropped as they always cancel with 
source terms in Tipi. 

2.5.1 the fermionic part 

Let us extract the fermion bilinears from eqs. (12. 3p . (12. 4p and (12.51) which are needed for 
our analysis in what follows. Rescaling the fermion fields so that their kinetic terms become 
canonical, we obtain 

-\ (g bb 9o b ; b n ^ b - \ (g bb g oh , b F°) v>V 

+ \ {^aag aa F°) \ a \ a ~ % - (F aag aa F°) A" A" 

-\ (g aa d a d a w) - \ {g aa dAW) rr. (2.16) 

Here the fermion fields ip a , ip a , A a , A a are to be integrated to make a part of the effective 
potential, while the gauge kinetic function J-* aa , the Kahler metric g aa and their derivatives are 
functions of the U(N) singlet c-number background scalar field </?°. The order parameters of 
supersymmetry F°, F°, and D° are taken as background fields as well. 
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From the lagrangian £p, the holomorphic part of the mass matrix is read off as 



-^FoaaF , v/^lHljy^Oaa^ 1 = | '"AA , 

■^y/g^iluLF^oaaD , g aa d a d a W + g aa g 0a ,aF° J ] 



We parametrize this matrix such that, in the case of F° = F° = 0, its form reduces to that of 
ref. [4.5J. The quantities having multiple indices such as J^oaa receive U(N) invariant expectation 
values: {J-o aa ) = (Fooo) e.t.c. See, for instance, [16]. We suppress the indices as we work with 
the unbroken U (N) phase in this paper. 

A S -*5«, (2.18) 

The two eigenvalues of the holomorphic mass matrix are written as 

A (±) = (tr.M)A {±) , (2.19) 

where 



A (±) = \\ 1 ± y (1 + iff + (l + I/) A^j . (2.20) 

These provide the masses for the two species of SU(N) fermions once the stationary values are 
determined. 

2.5.2 the bosonic part 

Likewise, we extract the bosonic quantum bilinears from eqs. (12. 3p . (12. 4p . and (12. 5ft . Let 

a = 5 o V + vV" 1 ^, (2-21) 



Ai = y/Qm7)-A%, (2.22) 
F a = ^g aa (ip)F a , (2.23) 
D a = v/(Im T) aa D a (2.24) 

where ip° are the background c- number field while (p a , A" F a and D a are the quantum scalar, 
vector and auxiliary fields respectively. 
We obtain 

_pa pa„v + papa + ™ 

4 2 



£g) = d^ a d^^* a - Lp a F aiiv + F a F a + -D a b a 



+F a ((^d a W) + (^AW)^) + F a ((V^W) + QT<9 a <9 a lW) . (2.25) 

We have also ignored — |(Re J r ) a b€ fil,fKT F^F^ as we eventually set ip a to be constant in our 
analysis and this term becomes a total derivative. 
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2.6 coupling to jV = 1 supergravity and super-Higgs mechanism 



If eq. (12.11) couples to N = 1 supergravity, the lagrangian is augmented to become the following 
one 1221: 



d 2 Q2£ 
1 



3 (VV - 8TZ)exp j-^[/i~($, $ f ) + r($, $ f , V)} | 



16^ 



+ /i.e. 



(2.26) 



The fermionic part of the lagrangian relevant to the super-Higgs mechanism is given by 



gd c T ab D a rX b 



2^2 



-^gd c *T* ab D^ c \ a - l -gD a i'^ a + \gD a ^\ a 



~g ah *D b *W*d a r cd \ c \ d - \g^D a Wd b ,T* cd \ c \ d 



(2.27) 



where e is the determinant of the vierbein and the covariant derivatives of several kinds are 
defined as follows, 



V,X = d,\ a + cu,X - gf abc A b \ c + -(d b KV^ b - d^KV^X 



b*\ \ a 



p 4 



+-gA b l ImF b \ a , 
D a W = d a W + (d a K)W, 
V a D b W = d a d b W + (d a d b K)W + 2(d a K)D b W - (d a K)(d b K)W. 



(2.28) 

(2.29) 
(2.30) 
(2.31) 



Now, we focus on the gravitino mass terms to discuss super-Higgs mechanism associated with 
eq. (E26D. 



gravitino mass 



y/2 [ V2 



i-^D a \ a + e K/2 D a W*i) a 



h.c. (2.32) 



The field redefinition of the gravitino 



V2 



V2 



(2.33) 



eliminates the mixing terms of the gravitino with the gaugino A and the adjoint fermion ip: 

e-^gravitino mass = -e^W^CJ^ + — -1— ^NG + h.C. (2.34) 



2W*e K / 2 

where the NG fermion ^ng absorbed in the massive gravitino is read 



^ NG = i4=A,A a + e K ' 2 D a *W*r. (2.35) 



The eq. ( 12. 34ft tells us that the gravitino mass is given by 

,(K)/2 



m 3/2 = e<*>'^. {2M) 



Requiring the cosmological constant to be almost vanishing 



0~(V) = ^(D") 2 + e K 



(2.37) 



the gravitino mass can be expressed in terms of the vev's of the auxiliary fields 



\(D a W)\*+$(D*)* 
\/3M| 



m 3/2 ~ e<*>/ 2 v - " - . (2.38) 



2.7 special cases 



As is mentioned in the introduction, the theory permits a sequence of interesting limiting cases. 
If we demand the Kahler function K to be special Kahler, K are expressible in terms of T as 

K = ImTr (2.39) 

and g a b = \va.jF a b etc. If we further demand such that the action possesses the rigid M = 2 
supersymmetry with one input function by choosing the superpotential to be a particular form, 
the tree vacua are shown to break M = 2 supersymmetry to M = 1 spontaneously pT5l4TTll23] 
0. We list the transformation laws for the doublet of fermions in this special case 



w-^u -^urn^- * 11 (240) 



where 



D a = -\g ab V b + ^=g ab [F bcd ^ d \ c + J^A C ) , 

F a = -V2Ng ab (e6° b + m.F 06 ) - % -g ab [T hcd \ c \ d - T hcd ^ d ) . (2.41) 



This superpotential consists of the terms referred to as the electric and magnetic Fayct-Iliopoulos terms. 
This M = 2 FI term is very special in the sense that, by the SU(2)r rigid rotation, it can be represented as a 
part of the superpotential. In this way, it avoids the difficulty (see, for instance, [24] for a recent discussion) of 
coupling the system to Af = 2 supcrgravity [231125] 
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2.8 connection with the previous work 

We here stop shortly to address the connection of ref. [I] with the previous work. Models of 
dynamical supersymmetry breaking with non-vanishing F- and D-terms have been previously 
proposed: they are, for instance, the 3-2 model [12] and the 4-1 model in [21]|^| In these models, 
supersymmetry is unbroken at the tree level and is broken by the non-vanishing vev of the F- 
term through instanton generated superpotentials. Non-vanishing vev of the D-term is also 
induced, but is much smaller than that of the F-term. 

In our mechanism, supersymmetry is unbroken at the tree level, and is broken in a self- 
consistent Hartree-Fock approximation of the NJL type that produces a non-vanishing vev for 
the D-term. A non-vanishing vev for the F-term is induced in our Hartree-Fock vacuum that 
shifts the tree vacuum and we explore the region of the parameter space in which F-term vev 
is treated perturbatively. 

We should mention that the way in which the two kinds of gauginos (or the gaugino and 
the adjoint matter fermion) receive masses is an extension of that proposed in [30]: the pure 
Dirac-type gaugino mass is generated in [30] u, while the mixed Majorana-Dirac type gaugino 
masse is generated in our case, the Major ana part being given by the second derivative of the 
superpotential. In [30] , the dynamical origin of non- vanishing D-term vev was not addressed. 

As for the application to dynamical chiral symmetry breaking, a supersymmetric NJL type 
model has been considered [I5HIE]- Chiral symmetry is not spontaneously broken in a super- 
symmetric case. Even in softly broken supersymmetric theories, the chiral symmetry broken 
phases are degenerate with the chirally symmetric ones. Thus, in supersymmetric theories, the 
phase with broken chiral symmetry is no longer the energetically preferred ground state. 



3 The effective potential in the Hartree-Fock approxi- 
mation 

The goal of this section is to determine the effective potential to the leading order in the Hartree- 
Fock approximation. We will regulate one-loop integral by the dimensional reduction [49] . We 



prepare a supersymmetric counterterm, setting the normalization condition. We make brief 
comments on regularization and subtraction schemes in the end of section 4. We also change 
the notation for expectation values in general from (...) to as our main thrust of this paper 
is the determination of the stationary values from the variational analysis. 



2 Application of these models to the mediation mechanism, see for example [27H29] , 
3 Attention has been paid to Dirac gaugino in many papers [251l!?!?ll5lTl44j . 
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3.1 the point of the approximation 

In the Hartree-Fock approximation, one begins with considering the situation where one-loop 
corrections in the original expansion in h become large and are comparable to the tree contri- 
bution. The optimal configuration of the effective potential to this order is found by matching 
the tree against one-loop, varying with respect to the auxiliary fields. In this section, we start 
the analysis of this kind for our effective potential. There are three constant background fields 
as arguments of the effective potential: (p = ip° (complex), U(N) invariant background scalar, 
D = D° (real) and F = F° (complex). The latter two are the order parameters of M = 1 
supersymmetry. 

We vary our effective potential with respect to all these constant fields and examine the 
stationary conditions. We also examine a second derivative at the stationary point along the 
constraints of the auxiliary fields to understand better the Hartree-Fock corrected mass of the 
scalar gluons. Let us denote our effective potential by V. It consists of three parts: 

V = V tmc + V c . t . + Vi_ioop. (3.1) 

The first term is the tree contributions, the second one is the counterterm and the last one is 
the one-loop contributions. After the elimination of the auxiliary fields, the effective potential 
is referred to as the scalar potential so as to be distinguished from the original V. 



3.2 the tree part 

To begin with, let us write down the tree part and find a parametrization by two complex and 
one real parameters. We also introduce simplifying notation goo(<p, <p) = g(<p, <p), (Im J r ((^)) o = 
Im F"((p),d W((p) = W'((f),g o,o = dg, etc. 

V tTee (D, F, F, <p, (p) = -gFF - U\mF")D 2 - FW' - FW> '. (3.2) 

As a warm up, let us determine the vacuum configuration by a set of stationary conditions at 
the tree level: 

Qytree 

-gg- =0, (3.3) 

Qytree 

= 0, as well as its complex conjugate, (3.4) 

oF 

Qytree 

—g — = 0, as well as its complex conjugate. (3.5) 
Eq. (13.31) determines the stationary value of D: 

D = = (3.6) 
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while from eq. (13. 4p . we obtain 

F^-g- x { V ,(p)W'{(p)=F^(p). 
Eq. (13.51) together with these two gives 



W'(<p*) = and therefore F*(ip, <p) = 0, 



(3.7) 



(3.8) 



as well as 



KELfo V) = V trCC (Vi <p,D. = 0,F= F.(<p, <p), F = F.(<p, (p)) = g-\<p, $)\W'(<p)\ 2 . (3.9) 

The negative coefficients of the RHS of eq. (13. 2 j) imply that both D and F profiles of the 
potential have a maximum for a given tp. These signs are, of course, the right signs for the 
stability of the scalar potential as is clear by completing the square. This is a trivial comment 
to make here but will become less trivial later. The mass of the scalar gluons at tree level |m s *| 2 
is read off from the second derivative at the stationary point: 



dipdip 



m s ((f,tp) =g^ l ((p,if)W"((p), 



(3.10) 

(3.11) 
(3.12) 



As we have already introduced in eq. (I2.18p . A and r are defined by 



A 



v^vVHlm^")- 1 ^" 



2 g-W + g~ x dgF 



D = r(tp,£,F,F)D. 



(3.13) 



Recall that we have suppressed the indices, invoking the U(N) invariance of the expectation 
values. Also 



/3 = 



where f% differs from / in eq. (12 . 1 8[) by 



g- v T"'F 



g- x W" + g- 1 dgF' 



(3.14) 



(g-'W" + g- l dgF)h = ( g^W + g~ l dgF - ~g~ x T"'F ) /. 



(3.15) 



We obtain 



f _j — ms a ]dg I f |2 



1 - 



(3.16) 
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While we will not make exploit in this paper, V tree can be written as a function of tp complex, 
|A| real, / 3 complex: 



V 



tree 



9 1 dg 

m s H : m s e 

g- l P" 



g j~ g~ -F 

(3.17) 



where m s ,e,g, J 7 (and their derivatives) are the functions listed above and undergo the varia- 
tions to be carried out in the subsequent subsections. We also see that the mass scales of the 
problem are set by m s *, the scalar gluon mass and g~ l T'l , the third prepotential derivative, 
(and g~ 1 dg), once the stationary value of the scalar is determined. 

3.3 treatment of UV infinity 

In the NJL theory there is only one coupling constant carrying dimension —1 and the 

dimensionless quantity is naturally formed by combining it with the relativistic cutoff, which is 
interpreted as the onset of UV physics. In the theory under our concern, UV physics is specified 
by the three input functions, K, J 7 , W and the UV scales and infinities reside in some of the 
coefficients. Our supersymmetric counterterm (HE] is 



Vet. = -^Im J d 2 9AW 0a W 0a = -^(lmA)D 2 . (3.18) 



It is a counterterm associated with ImJ 7 ". We set up a renormalization condition 

1 d 2 V 

= 2c, (3.19) 

D=0,(p=tp*,ip=(pt, 



N 2 (dD) 2 



and relate (or transmute) the original infinity of the dimensional reduction scheme with that 
of ImJ 7 ". We have indicated that this condition is set up at D = and the stationary point of 
the scalar which we will determine. We stress again that the entire scheme is supersymmetric. 

3.4 the one-loop part 

The entire contribution of all particles in the theory to %• (the 1PI to one-loop) = iY\_\ oov is 
easy to compute, knowing ( I2.19p . f)2.20p and It is given by 

■r ( fj' W f d ' k i AiAi +> l 2 -* 2 -fe)(iAi~'l 2 -fc 2 -^A 
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In the unbroken U(N) phase, it is legitimate to replace by iV 2 and drop the index a as we 
have said beforej^ We obtain 

, , 1 



loop : 



1 } (fd*x 

-iV 2 |tr.M| 4 

iV 2 |tr.M| 4 J. 



r 



1— loop 

dH» 



(3.21) 



In 



( {\\^\ 2 -l 2 -ie){\\^\ 2 -l 2 -iej 



v (i^&r-p-fe)H 



— ie\ 



(3.22) 



Note that |m s | 2 , whose stationary value give the tree mass squared of the scalar gluon, differ 
from |tr.M| 2 : 

2 



\XxM\ A 



m s --(g- 1 F")F+(g- 1 dg)F 



(3.23) 



To evaluate the integral in d-dimensions, we just quote 

dH" 



I(x 



I{x 2 ) - 1(0) 



(2tt)H 



log(x —I —is) 



1 



32vr 



[A{e^)(x 2 ) 2 -(x 2 ) 2 



log(x 2 )] 



where 



. / N. l l 



e = 2 



d 



(3.24) 
(3.25) 

(3.26) 



We obtain 



loop 



N 2 \trM\ 
32vr 2 



A(e,i) |A (+) | 4 + |A(")| 4 



m 



-|A( + )| 4 log|AW| 2 -|A(-)| 4 log|A(-f + 
This again depends upon A, / and <p. 



tiM 



tiM 



log 



tiM 



(3.27) 



4 In those cases where the U(N) is broken to the product group flLi we nee d n °t only replace 

• • • by • • -q, but also must treat the Af = 1 multiplct of the broken generators that receives the mass by 

a a 

the Higgs mechanism [16] . 
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4 Stationary conditions and gap equation 



4.1 variational analyses 



Now we turn to our variational problem. It is stated as in the tree case as 

3D ' 
dV 

— — = and its complex conjugate, 
oF 



8V_ 

dip 



and its complex conjugate. 



(4.1) 
(4.2) 
(4.3) 



We will regard the solution to be obtained by considering eqs. (14. ip and H4.3[) first and solving 
D and tp for F and F: 



D = D*(F,F), <p = <p*(F,F), <p = <p*(F,F). 

Eq. is then 

dV(D = D*(F, F),tp = tp*(F,F),<p = <P*(F, F), F, F) 



OF 







(4.4) 



(4.5) 



D,ip,<p,F fixed 



and its complex conjugate. These will determine F = F*,F = F*. 

In this paper, we are going to work in the region where the strength is small and can 
be treated perturbatively. This means that, in the leading order, the problem posed by eq. 
(14.11) and eq. (14.31) becomes 

dV(D,<p,p,F = 0,F = 0) 



dD 



0. 



dV(D, cp, (p, F = 0, F = 0) dV(D, <p,<p, F = 0, F = 0) 







(4.6) 
(4.7) 



dip dip 

and this problem does not involve the tree potential eq. (13. 2 p except the last D 2 term, as F, 
and F are set zero. Eq. (14. 6 p is nothing but the gap equation given in [HE], while eq. (14. 7p is 
the stationary conditions for the scalar. This is the variational problem which is analyzed in 
this paper. A set of stationary values (D*, <^*) is determined as the solution. 



4.2 the analysis in the region « 

Let us first determine V(D, ip, (p, F = 0, F = 0) explicitly. We need to solve the normalization 
condition. 



2cN 2 



d 2 V 
(dD) 2 



-(ImF") - (ImA) + A^ 2 |trA4|^ 



D=0,« 



d 2 J 

(dDf 



D=0 
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where J has been introduced in eg. (13.221) . At F, F — > 0, 



where 



A — > Ao = To((f, y?)-D, where tq 



V2 y/g-\\mP')-^r 
2 g- x W" 



x(+) . \(+)\H 



X A 2 A (+) A ( ~) 



1 + Ag, 



(4.9) 
(4.10) 

(4.11) 



tiM 



-> 1, 



J J 



327T 2 



A(e )7 ) 



|Ai +) | 4 log|A^T-|An 4 log|A^ 



-A 



,(+)|2 



H|4l 



1 + AnA/l + An— 1 



,(")|2 



(4.12) 



(4.13) 



essentially reducing the situation to that of refs. [IJ|5]. 

Note, however, that r and A (or tq, Aq) are complex in general except those special cases 
which include the case of the rigid M = 2 supersymmetry partially broken to M = 1 at the 
tree vacua. For |A | <C 1, 



Jo 



32tt 2 



A(e n )l(Al + At)-±(Al + At) + O(\A \ 



We solve the normalization condition for the number A to obtain 



.4: 



+ 



327T 2 



2 ' \m a >\*{rl+rl) 



2c 



ImJ 7 " ImA 



N 2 



We obtain 



N 2 



(4.14) 



(4.15) 



V o = V(D,<p,<p,F = 1 F = 0) 
= -hmF"D 2 - ^(ImA)D 2 



+ 



iV 2 |m, 



327T 2 



A(c, A, <p*,(p* 



-iAi + 'riogiA< +) r-iA<-'riogiA<-)| 2 



1 + 5 A » 



After some calculation, this is found to be expressible as 

2 



V 



( 1 



N 2 \m s 



64^ + S -^'^ 



A + A 



l + A 2 Vi + A 2 -l 



(4.16) 



A -|A | 4 + /(A ,A C 



t (-) |4 



32 ^(|Ai +) | 4 log|Aj +) | 2 + |An 4 log|A 

16 



32tt 2 V 8 

H|2 



(4.17) 



where 



14 I r L+rl 



%> <P) = 7i 



/(A ,Ao) 



Im^y _,_ ImA 

Af 2 A^ 2 

MilfoJl™ 14 
2 



ImJ r/ 7A r2 +ImA/Af 2 
Im^'/Af 2 +ImA/Af 2 



\m s \ 4 



, r o+ r o ' 



|m,.|' 



A 2 



A 2 



lAnl 2 - 1 



(4.1* 



(4.19) 



(4.20) 



Note that 



5* = 5(y?*,y?*) ^ 0, 



(4.21) 



and 



|/(A , A )| < const for |A | > 1. 



(4.22) 



If r (and A ) is real, this is rewritten as 



V 



iV 2 |mJ 4 



64vr 



5 I An 



327T 2 



^-A« 4 logAW 2 



(4.23) 



where c' 



is the rescaled number, and 



1 / lmT'J | ImA 

1 / Af2 i- ^2 

2 y ?"o*l m s*| 4 



ImF" /N 2 +ImA/N 2 
ImJ'i'/A^+ImA/Af 2 



(4.24) 



and m a ((p,(p) = g~ 1 W" are the functions of </?, Clearly, there are two scales in our current 
problem |rn*| -1 / 2 and |m s *|, which are controlled by the second superpotential derivative and 
the third prepotential derivative at the stationary value y?*. 
Let us turn to the gap equation 



dVo 
dD 



0. 



(4.25) 



tp,tp 
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For eq. (14.17[) . scaling out |r | 2 , we obtain 



64vr 2 



+ c-8 (1 + cos 20) + 



.4 



32vr 2 2 



;|A | 2 - (1 -COS2(0-0')) 



1 



2 log I A 



(+)|2 



327T 2 

- (21og|A - ) | 2 + 1 



1 / e 2W X 



2i6 \ (+) 



2i6 \ (+) 







+ 



e~ 2W X 



o 



2 VVi + Ag v^Ta 2 



(+)|2 




1 / e 2M X 



2i6 \ (-) 







2i9\{~) 



e~™X 







2 VaA+^I v/l + A 



H|2 



where 



A = |Ao|e tfl , r =|r |e lfl , tan 20' 



|A o | 2 sin20 



(4.26) 



(4.27) 



1 + |A o | 2 cos20' 

Note that |1 - cos 2(0 - 0')| ->■ in the region ~ or |A | > 1. 

On the other hand, for eq. (I4.23P with Ao being real, N 2 \m s \ 4 is scaled out and it is simply 
given by the An derivative: 



= A 
= A 
+ 



c' + 



64tt 2 

J-) 

64vr 2 J 



4a2 + 1 d 



32tt 2 1 2 ' A dA 



-AriogAr-AS- )4 iogA(-) 2 



32tt 2 



1 



'ITA 2 



A n +)3 (2 log A 



(+)2 



Al ( 2 log A 



(-)2 



(4.28) 



which is our original gap equation^ In both cases, the solutions are given by the extremum 
of the potential Vq(D, tp, <p) in its D profile. We stress again that the D profile is not a direct 
stability criterion of the vacua, which is to be discussed with regard to the scalar potential 
V o (D m ((p,0),tp,0). 

and its complex conjugate. For eq. (I4.17p . we obtain 



We next examine ^ 



D,tp 



d 

2—— (In m s 

dp y 1 



V 



N 2 \m s 



d6_] (A + A 

dp 



<9 In rn 


H 


din rn 




d ( V \ 


dp 


<p 


dp> 


p. 


dD \N 2 \m s \ 4 J 



-DV 



Vn 



N 2 \m s 



D 



and its complex conjugate 



where 



d 



d 



'dAr 



r 



dAr 



(4.29) 



(4.30) 



5 We have introduced 5(ip, ip) such that its stationary value 5(<p*, tp*) =0, which can therefore be ignored in 
analyzing eq. (|4.28l) . 
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The second term of the RHS of eq. (14.291) is proportional to the gap equation eq. (14.26)) . As 
for the third term, after some calculation, we obtain 



-V 



Vo 



JV 2 |m B | 4 



f )|A||ro| 



64vr 



+ c - 5 sin 29 



32tt 



:Asm2(8-6') 



1 UBin(2^-^ + sin2(0 _ 0/))|A ( +) w 21og|A ( +)|2 + i 



32vr 2 2 V |1 + Ag|V: 
1 1 /sin(20 - 9' 



32tt 2 2 
:C(0,|A o |) 



(iTTAlW - sin2 ^ " e ')) l^l 2 ( 21 °s l^l 2 + 1 



(4.31] 



In the RHS of eqs. (I4.29P and (14.311) . we have regarded A , A , (p and (p as independent variables. 
For eq. (14.23)) . with A real, we obtain 



29(ln \m s 



Vo 



<9A n d 



N 2 \m q 



f d5 \ A 2 



v 



dip 8A \N 2 \m s \ A 



(4.32) 



and its complex conjugate. Here in the last term of the RHS, we have regarded Ao, <p, <p as 
independent variables. 




0.0 0.2 0.4 0.6 



Figure 1: The schematic picture of the intersection of the two curves which represent the 
solution to the gap equation (the red one) and the p flat condition (the blue one). The 
horizontal axis is denoted by p/M and the vertical one by Ao- The values at the stationary 
point (A *,y* = <p*) are read off from the intersection point. 

Finally the stationary values (D*, </?*, </?*) are determined by eqs. ( 14.26)1 and ( 14.29)1 or by 
eqs. (I4.28P and (I4.32p . Let us discuss the latter case first. As the second term of the RHS in 
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eq. (I4.32p is nothing but the gap equation eq. (14.281) . eq. (I4.32[) can be safely replaced by 



V dip a 2 



as 

K (4-33) 



N 2 \m s \ 4 2<9(ln|m s | 2 

<p being real. The solution to eq. ( 14 . 3 3 1) in the A profile is determined as the point of 
intersection of the potential with the quadratic term having tp — (p dependent coefficients. 
Actually, it is a real curve in the full (Ao, tp — (p) plane. Likewise, the solution to the gap 
equation eq. (I4.28p . the condition of A extremum of the potential, provides us with another 
real curve in the (A ,y? = <p) plane. The values (A *,y?* = tp*) are the intersection of these 
two. The schematic figure of the intersection is displayed in Figure [U By tuning our original 
input functions, it is possible to arrange such intersection. Conversely, as an inverse problem, 
for given A * and the height of the A profile, one can always find the values of the coefficients 
in eq. (14.281) and the coefficient function in eq. (I4.33P that accomplish this. Dynamical 
supersymmetry breaking has been realized. 

As for the former case, as in the latter case, we can safely replace eq. (14.291) by 



95 \ - ■ —2 (M) |A | 2 



Vq [dtp) /A + A \ , \8 V 



-C. (4.34) 



iV 2 |m s | 4 2<9(ln|m s | 2 ) \ 2 J 2d(\n\m 

The values (A *, A *, tp*, tp*) can be determined by the intersection of eq. (I4.26P and eq. (I4.34p . 
We will not carry out the (numerical) analysis for this case further in this paper. 

4.3 determination of 

Let us now turn to the analysis of the remaining equation of our variational problem, eq. (14. 2p . 
In our current treatment, 



F = --W' + -7^1-loop « --W + -iv^op 

9 gdF g gdF 



(4.35) 



F=0 

As the stationary values (D*, are already determined, this equation and its complex 

conjugate determine F* and F: 



-W (<p*, <p*) + —=;Vi-i oop (D*,ip*, <p*, F, F) 



g{<P*,tp*) V ' 9F 



F=F=0 



(4.36) 



Note that, knowing Vi_i o P explicitly in eq. (I3.27p . the RHS can be evaluated. We can check 
the consistency of our treatment through f 3 in eq. (I3.14p by I/3I 1. 
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4.4 numerical study of the gap equation 



In this subsection, we study some numerical solutions to the gap equation eq. ( 14.280 and the 
stationary condition for p eq. ( 14.330 in the real Aq case. The equations we should study are 








Vn 



A 2 



,(+)3 



dip 



N 2 \m K 



29(ln |m s „,| 2 



21ogA« 2 



<P*,<fi* A 2 



1 - A 



(")3 



21ogA^ )2 + l 



A =A * • 



(4.37) 



(4.38) 



where we note that 8(tp, (p) in the gap equation ( 14.280 vanishes at the stationary point in the 
real Ao case. By using eqs. ( 14.230 and ( 14.240 . the second condition can be rewritten after 
dividing by as 



c' + 



64tt 2 
1 



+ 



A 



32vr 2 
Im^' + A) 



A 2 



A 2 



A 



(+)|4 



log I A 



(+)|2 




I A 



(")|4 



log|A^| 2 



=A , 



4A 2 <91n 



o* 



din Im(J r,, + A) | 



<p*,v>* 



d(r 


m s \ 


2^2 


* 




Wis* 


4 ) 2 



(4.39) 



The nontrivial solution A * ^ to the gap equation (14.370 is found by some region of the 
parameters d and A, which was already done in [I]. This solution fixes the LHS of eq. ( 14.390 
and y?* is determined by solving eq. ( 14.390 in principle. In order to find ip* explicitly, the form of 
the prepotential J 7 and that of the superpotential W must be specified. Here, we take a simple 
prepotential and a superpotential of the following type (with some abuse of the notation): 

,3 



-r- C n 1 

T= tro? 2 + 

2N y 3\MN 



m 2 d 
W = — troj + — — 
N * 3\N 



tip d 



W = \cp 2 + -L^, 
d 



™> 2 <p + t^ 3 , 



(4.40) 
(4.41) 



where c, d are dimensionless constants while m, M carry dimensions. In particular, M is a 
cutoff scale of the theory. This prepotential is minimal for DDSB. As for the superpotential, 
at least two terms are required to be supersymmetric at tree level. We can take a quadratic 
term (p 2 instead of the cubic one, but in that case, RHS of eq. ( 14.390 becomes singular because 
of <9 In ImJ 2 = 0. 



Substituting these T and W into eq. (14.390 . we obtain 



c + 



64vr 2 



32vr 2 



A A2 
—A 



o* 



A 2 



,(+)|4 



log I A 



(+)|2 



IA 



(")|4 



log I A 



H|2 



Im(c + A)(Im c) 



A 2 



(d<p*/M) 



2 ■ 



(4.42) 
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where we utilized the fact that 1/M, d, tp* are real and c is pure imaginary, which are necessary 
for A = A . If we take the coefficients c — i,d — 1 for further simplification, we can easily 
obtain a solution by tuning N and ImA. We note < <p*/M < 1 for our effective theory to be 
valid. In our analysis carried out in this paper, we consider the region where the magnitude of 
the F-term is smaller compared to that of the D-term. Therefore, we need to check whether 
our solutions satisfy this property consistently. Let us consider the ratio of the auxiliary fields: 

-0 -1 W(£,) + g" 1 JfVi-ioop^*, f*, F, f )\ f= f=o 



D 



+i 



+ 



N 2 



/ m 



tp, 



x ~ 2 



Ao* / ro* 
2 \M) 



^-A 2 

1287T 2 °* 



32tp 



: (|A + J 4 log|A + J 2 +|A -| 4 log|A 



o* I 



327r\/l 



A 2 



(A + J 3 (log | A 



, 2 + 
0*1 



(A 



0*/ 



log|AoJ 2 + 



(4.43) 



where the form of the prepotential and that of the superpotential in eq. (I4.40p and eq. ( I4.66P 
are assumed and we have put c = i, d = 1 in the second equality. 

Now, the numerical solutions to the gap equation and the stationary condition for tp are 



listed in Table [TJ In these examples, we have taken some values of 



N 



and m just for 



Im(i+A) 

an illustration and the ratio \F*/D*\ and 1/3*1 are evaluated. We can find that the F-term is 
smaller than the D-term in some of these examples. 



4.5 second variation of the potential and the mass of the scalar 
gluons 

We now turn to the question of the second variations of the scalar potential 

Scalar = V(D = D.(<p, tp), F = F*(<p, tp) w 0, F = F,{<p, tp) « 0, tp, tp) {AAA) 

at the stationary point (D*(<£>*, tp*), 0, 0, </?*, </?*). It is convenient to separate V(D,F,F,ip,tp) 
into two parts: 



Here 



V = V + V . (4.45) 



V(F, F, if, tp) « —gFF - FW' - FW' + {dpVx-^F + {dpV^^F 

+\(d 2 F V l ^ oop )*F 2 + ^{d 2 F V^ loop )*F 2 + (d F d F V^ loop )*FF, (4.46) 
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J 1 1 


i/(4 • 32vr 2 ) 


A * 


<P*/M (-^) 


\F*/D m \ 


1/3*1 


C + 64tt 2 




0.002 


0.0001 


0.477 


0.707 (10000) 


2.621 (m = M) 


1.77 


0.002 


0.0001 


0.477 


0.707 (10000) 


0.524 (m < M) 


0.35 


0.002 


0.0001 


0.477 


0.707 (10000) 


0.860 (m = 0.4M) 


0.58 


0.003 


0.001 


1.3623 


0.8639 (2000) 


0.825 (m = M) 


>1 


0.003 


0.001 


1.3623 


0.8639 (2000) 


0.224 (m < M) 


0.43 


0.003 


0.001 


1.3623 


0.5464 (5000) 


1.092 (m = M) 


>1 


0.003 


0.001 


1.3623 


0.5464 (5000) 


0.142 (m < M) 


0.27 


0.003 


0.001 


1.3623 


0.5464 (5000) 


0.911 (m = 0.9M) 


1.76 


0.003 


0.001 


1.3623 


0.3863 (10000) 


1.444 (m = M) 


>1 


0.003 


0.001 


1.3623 


0.3863 (10000) 


0.100 (m < M) 


0.19 


0.003 


0.001 


1.3623 


0.3863 (10000) 


0.960 (m = 0.8M) 


1.85 



Table 1: Samples of numerical solutions for the gap equation and the stationary condition for 
if. The ratio \F*/D*\ and 1/3*1 are also evaluated for consistency check. 

and 

V (D, <p, <p) = V(D, ip,if,F = 0,F = 0). (4.47) 

In eq. ( 14.4611 . we have extracted the F,F dependence of Vl-1oo P (eq. (13.2711 ) as series and * 
indicates that they are evaluated at (£>*, f*,(p*, 0, 0) after the derivatives are taken. Eq. ( I4.47[l 
has been computed in eq. ( I4.16P and eq. ( I4.23p . We will compute the second partial derivatives 
and the second variations of Kcaiar, using the formula in the appendix. 
For V, y L = {F,F),y R = (ip,(p), 

( d 2 V, ddv\ 

M -H^;^j^°' (4 - 48) 

M ^fdd F V,dd F V\ ^ f-W" + (dd F V 1 ^ loop ), _(dd F Vi-y oop ) \ 

RU \dd F V,ddpV) ~\ {Bd F V^ loop ), -W" +{dd F V 1 -i oop )) ^ 
M L R. = M t RL ., (4.50) 

M / d F v, d F d F v\ / (sivs-ioop), -^ + (^y 1 _ loop )\ 

LL * \d F d F V, d 2 F V )~ \-g + (d F d F V^ loop ), (a|T4-ioo P ) 

Here we have denoted by * that the derivatives are evaluated at the stationary point. 
We obtain, after some computation, 

<m - W R M R U-MEljM LR Jy R = \tjf R (^f^/A Sy R . (4.52) 
2 2 KM^M^J 
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Here 

= g[G 2 l _\ C?) ( G (^l 2 + I^H + CAB + CAB ) > ( 4 - 53 ) 
Mw = g{G 2 -\C\ 2 ) ^ GAB + + ' (454) 

g ^ t _ ^Vl-locp A ^ w „ _ ddp y B = -ddpV^. 

9 9 

(4.55) 

Here in the last line of eq. fl4.52p . we have changed the real quadratic form into the complex one. 
We see that in the region |(<9i7<9^V)o|*, |(<9f^) |*, <C g*, the matrix J\4* is well approximated by 

1 / \A\ 2 + \B\ 2 , 2AB \ 
M* ~ - ( - - . .„ , „ • 4.56 
g\ 2AB, \A\ 2 + \B\ 2 )^ 

The two eigenvalues are 

~{\A\ ± \B\)l = - {\W" - (dd F V^ loop )\ ± \(dd P V^ oop )\)l , (4.57) 
9 9 



respectively, ensuring the positivity of ( 14 .5^ 

For V , y L = D, y R = (ip, ip), 



f d 2 V , ddV \ I dd D V \ 

Mrr = ' _ , M RU = , 4.58 

\ddVo,&Vo)j \ddnVoJJ 1 J 

M LR , = M* RL , M LL , = d 2 D V *. (4.59) 

We know that the D profile of Vq(D, tp, (p) near the stationary point is convex to the top and 
we fit this by 

V = V h (<p, ip) - (D - D m (<p, p)f + 0{{D- D*(<p, p)f). (4.60) 

Here a is a positive real function of <p, <p and Vh(<p, <p) = Vb(-D*(v?, ip), <p, ip). One can check 

- M RL Ml l L M LR * = a* I dB * } ((dD*), (3D*))*, (4.61) 



while 



3D* 



MRR * = { 8 dVh 5 2 v h - a [\6D.nBD.r ] (462) 
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and 



<5V * = -5y R (M RR * - M RL *M£^M LR *)Jy R 




ddV h d 2 V h 



R = 5y R M h *y R . 



(4.63) 



d 2 V h ddV h 

■ * 

The entire contribution of the second variation 5 2 K = <5 2 V* + <5 2 Vo* to the leading order in 
the Hartree-Fock approximation is given by eqs.f l4.52j) . (14. 63ft . The mass of the scalar gluons 
squared is obtained by multiplying the combined mass matrix by g~ x : 

g^iMt + Mk*), (4.64) 

generalizing the tree formula. In practice, we just need a well-approximated formula valid in the 
region we work with and one can invoke the U(l) invariance to ensure that the two eigenvalues 
of the complex scalar gluons are degenerate. Let us, therefore, use the expression 

|2 



-\w" 

9 



^-ziV 2 
M 

1 



(dd F Vi 



1— loop J 



ddV } 



h* 



M 



A(e,i) 
32tt 2 



(Ai) 4 + (A ") 4 + 



(p*/M 



-A 



o* 



-A 4 



32tt 2 



{(Ai) 4 log(A+) 2 + (A -) 4 log(AoJ 2 



3(1 + % 



(A - r log(A 



+2N 



M) 
Im(i + A) 



647P 



A 2 



A 




32tt 2 \ 8 



t ~ (AoJ 4 log(A c J;; 



32tt 2 cp*/M 
(AoJ 4 log(A -^ 2 



M 



jv wmj^ (465) 

to check the local stability of the potential and the mass. The above expression is obtained for 
our simple example of T and W 



1 



w 



(4.66) 



2 r ' 3!M 

and the real case A = A is applied. Using the numerical analyses carried out in the last 
subsection, we have made a list of data on eq. (I4.65p . Except for the last case in the table 
[21 the scalar gluon masses squared are found to be positive for any N, which implies that our 
stationary points are locally stable. Even in the last case, the stability is ensured for small N. 
In these data, we have checked that the inequalities | (dpdp V\ _i 00 p) |*, |(<9|<Vi_ioop)|* *C g* are 
in fact satisfied. As a summary of our understanding, a schematic figure is drawn in Fig. (2j 
which illustrates the local stability of the scalar potential at the vacuum of dynamically broken 
supersymmetry in comparison with the well-known NJL potential. 
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i/(4 • 32vr 2 ) 


A * 




scalar gluon mass 


C + 64tt 2 




0.002 


0.0001 


0.477 


0.707 (10000) 


0.4998 + 0.0056 N 2 + 8.607 x lO^N 4 


0.003 
0.003 
0.003 


0.001 
0.001 
0.001 


1.3623 
1.3623 
1.3623 


0.8639 (2000) 
0.5464 (5000) 
0.3863 (10000) 


0.7463 + 0.0106 N 2 + 2.653 x 10~ 4 iV 4 
0.2986 + 0.0008 N 2 + 4.694 x 10" 5 iV 4 
0.1492 - 0.0024 N 2 +7.235 x 10" 5 iV 4 



Table 2: Samples of numerical values for the scalar gluon masses. 




D 



Figure 2: Comparison of V^ ca i ar around the stationary value (D*,ip*) with Vnjl- 

4.6 summary and choice of regular izat ion and subtraction scheme 

In this paper, we have considered the theory specified by the general M = 1 supersymmetric 
lagrangian eq. (12. ip . have regularized the theory by the supersymmetric dimensional regulariza- 
tion (dimensional reduction) and have subtracted the part of the 1 / (e) poles of the regularized 
one-loop effective action in eq. ( 13. 27ft by the supersymmetric subtraction scheme defined by the 
condition eq. ( I3.19p . The upshot is an effective potential eqs. ( I4.17p . (I4.23P as a function of the 
background constant scalar and the order parameter D of supersymmetry, with another order 
parameter F of supersymmetry being induced and treated perturbatively. Supersymmetry is 
dynamically broken as is represented by the non-vanishing value of the order parameters at 
the stationary point. The original infinity is transmuted into the infinite constant A which is 
the coefficient of the counterterm and the effective potential has been recast to describe the 
behavior of the theory well below the UV cutoff residing in the prepotential function. As the 
theory is perturbatively nonrenormalizable, A is still present in our final expressions of the 
effective potential and we regard it to take a large value. 

We now make brief comments on other regularizations and subtraction schemes which we 
did not employ in this paper. The relativistic momentum cutoff is a natural choice of the 
NJL theory as we mentioned earlier but regularizing the integral eq. (I3.20p by the momentum 
cuoff leads us to a rather unwieldy expression. See ref . [5] . Unlike supersymmetric dimensional 
reduction [49], the momentum cutoff perse, while preserving the equality between the bose 
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and fermi degrees of freedom, does not have a firm basis on the regularized action which the 
supersymmetry algebra acts on. Moreover, as is clear from (Al) of ref. [5], the result violates 
the positivity of the effective potential in the vicinity of the origin in the A profile. This 
violation is a necessity in the broken chiral symmetry of the NJL theory but here it contradicts 
with the positive semi-definiteness of energy that the rigid supersymmetric theory possesses. 
Turning to the choice of the subtraction scheme, one might also like to apply the "(modified) 
minimal subtraction scheme" in our one- loop integral eq. ( 13.271) . While we do not know how 
to justify this prescription here, the subsequent analyses proceed almost in the same way and 
the main features of the equations obtained from our variational analyses and the conclusions 
are unchanged. 
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A NJL effective action 



In this appendix, we briefly recall a few aspects of the SU (N) Nambu-Jona-Lasinio model. 

£ NJL = Wfi& + ^ [(# *) 2 + (m^) 2 } (A.1) 
The equivalent lagrangian is 



£ = £ 



1 1 



NJL 



2X/N 



2A/iV 2X/N 



(A.2) 



The 1PI vertex function (or the effective action) Txpifcr, n] to one-loop (or leadig order in 1/N) 
reads 



«Tipi[<7, 7r] 

The gap equation is 



i 1 



2 A/iV 



d 4 x (a 2 + 7T 2 ) +N In Det (i# - a(x) - i^ir(x)) . (A.3) 







iST 



ipi 



5cr(x) 



.1 AT f d 4 k 1 

— —l ; (Tn — JM I —ri-T- 

CT (,)=H-o,^)=o X/N J (2tt) 4 ^-a 



1 f ^ 

■ cr I -Z-^tt - 4iV / — — X 



X/N 



(2vr) 4 fc 2 - a 2 



(A.4) 



B formula for the second variation 



In this appendix, we recall the formula for the second variation of a multivariable function 
subject to a set of stationary constraints. Let V be the function of two sets of variables: 
{{y\ • ■ ■ , u n ^}} = V L , {{y n{L)+1 , ■ ■ • , y»( L )+*( a )}} = V R . Namely, 



under 



V = V(y\--.,y^ L \y^ + \---,y^ + ^) 



dV 

— = ,i = 1,... , n (L). 



Let the second variation of V be 



d 2 V 



s 2 v = - y 

2 dy l dyi 

y\yj£V L UV R 



Sy % 6y j 



(B.l) 



(B.2) 



(B.3) 



but Sy l e T>l are not independent variations. 
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It is convenient to introduce a new vector notation: 



y L =(y\... ,y^)\ y R = (y 



n(L)+l . _ . n(L)+n(R)\t 



V 



] )\e.t.c 



(B.4) 



Define 



d 2 V 
dyidyi 



y i E V x , y j e V x , X, X' are either L or R. (B.5) 



Eq. flB~3]) reads 



1 1 
$ 2 V = -(5y R , M RR 5y R ) + (Sy R , M RL 5y L ) + -(5y L , M LL Sy L 



(B.6) 



Varying (1B.2|) with respect to yi and y R , we obtain 



M LL 5y L + M LR 5y R = 0. 



(B.7) 



Hence 



II = -M7}M 



LL m LR oy R . 



(B.8) 



Substituting this into (1B.6j) . we obtain 



S 2 V = l - (5y R , (M RR - M RL MllM LR ) 6y R ) . 



(B.9) 



The generic scalar mass matrix in the text can be read off M rr —M R lM l ^Mi R at the stationary 
value. 
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